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Abstract 

A new method is proposed for extracting limits on spin-dependent WIMP-nuclcon inter- 
action cross sections from direct detection dark matter experiments. The new method has 
the advantage that the limits on individual WIMP-proton and WIMP-neutron cross sec- 
tions for a given WIMP mass can be combined in a simple way to give a model-independent 
limit on the properties of WIMPs scattering from both protons and neutrons in the target 
nucleus. Extension of the technique to the case of a target material consisting of several 
different species of nuclei is discussed. 
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1 Introduction 



Weakly Interacting Massive Particles (WIMPs) are believed to be the most plausible candidate 
for dark matter (DM) in the Universe. WIMPs are predicted to exist in many extensions of 
the Standard Model of particle physics. Most of the well- motivated WIMP candidates are 
Majorana (i.e., X = X) fermions. This in particular is often the case in models based on 
supersymmetry (SUSY). Perhaps the most popular WIMP candidate is the lightest neutralino, 
a superposition of the SUSY partners of the electroweak gauge bosons (gauginos) and Higgs 
particles (higgsinos). Other plausible WIMP candidates for cold DM (CDM), the axino and 
gravitino (fermionic SUSY partners of the axion and graviton, respectively) are also Majorana 
particles. 

In the case of non-relativistic Majorana WIMPs, their predicted elastic scattering couplings 
to atomic nuclei are effectively of two types [jjj. In scalar, or spin-independent (SI) interactions 
the WIMP coupling is proportional to the mass of the nucleus. In the axial, or spin-dependent 
(SD) case, the coupling is proportional to the spin of the nucleus. 

Direct searches for galactic halo WIMPs through their elastic scattering off target nuclei 
are currently being carried out by several groups. In the absence of a positive signal it is the 
aim of these experiments to set limits on the properties of WIMP dark matter independently 
of its precise composition. This is accomplished by setting limits on the cross sections for 
SI and SD interactions between WIMPs and target nuclei as functions of WIMP mass. To 
enable comparison with results from other experiments (which may use different target nuclei) 
one often translates these limits to bounds on the WIMP-proton cross section. The current 
procedure used for the conversion is relatively straightforward in the case of SI interactions 
but becomes problematic in the case of SD interactions. 

The reason for this is that limits on spin-dependent WIMP-proton scattering cross sections 
contain considerable dependence on a particular WIMP composition {e.g., gaugino-like versus 
higgsino-like neutralino WIMP). Since the spin of target nuclei is carried both by constituent 
protons and neutrons, when converting to a WIMP-proton cross section a value for the ratio of 
the WIMP-proton and WIMP-neutron cross sections must be assumed. But in many cases this 
ratio can vary significantly depending on the assumed type of WIMP. In the particular case of 
a predominantly gaugino neutralino WIMP this ratio varies by several orders of magnitude § . 
As a result, current experimental limits on the WIMP-proton cross section for SD interactions 
are fraught with potentially significant WIMP-type dependence. 

In this Letter, we introduce an alternative method for deriving limits on spin-dependent 
WIMP-nucleon interactions. This new procedure allows one to extract experimental limits 
in a (Majorana) WIMP-model independent way. First, in Sec. 2, we briefly summarise the 
current practice. The new method is presented in Sec. 3 and its features are discussed in 
Sec. 4. 

2 The Current Procedure 

The current procedure for calculating spin-dependent WIMP-nucleon cross section limits from 
experimental data can be summarised as follows (following Refs. [||, H|). Assume a detector 
consists of some species of nucleus, = N . Using the convention of Ref. Q (see Eqn. (7.13)), 



1 



the total WIMP-nucleus cross section a a can be written as 



a A = AGI&Ca, (1) 

where the WIMP-target reduced mass \ia is given by m x mA/(m x + tua) for WIMP mass m x 
and target nucleus mass tua- The "enhancement factor" Ca is given in Eqns. (7.14) and (7.35) 
of Ref. for SD and SI interactions, respectively, and will also be given explicitly below. 

As emphasized in Ref. []|], a a (called there <ro) is, strictly speaking, not the total cross sec- 
tion. It is, by definition, the "standard" total WIMP-nucleus cross section at zero momentum 
transfer. However, it is the quantity that is conventionally used by experimental groups for 
setting limits and we will continue calling it here the total WIMP-nucleus cross section. 

In the first step, the data is used to calculate limits a A m on the SI and SD WIMP-nucleus 
cross sections. (If the target is made of several species of nuclei (e.g., Na and I), the procedure 
is performed for each species separately and next the limits are combined together.) In each 
case it is assumed that only the given type of interaction (SD or SI) dominates the total cross 
section. 

The WIMP-target cross section limit a\ m obtained for target A can be expressed in terms 
of limits on WIMP-nucleon (i.e. free proton or neutron) cross sections <7p im ^ and ctvJ A ' 

2 2 
lim(A) lim ^P 1 ,Jim(A) lim - 1 / \ 

a * =(Ja ACa1Cp n aA A~cIiCn (2) 

where /i P) „ is defined by setting A = p (A = n) in the expression for the reduced mass ha 
above, and similarly for C p ^ n (see also below). This conversion is made conventionally to the 
WIMP-proton cross section limit using the former expression. 

The purpose of this conversion is twofold. It allows one to compare limits derived by 
different experiments which use different target materials. Second, theoretical calculations 
in specific (e.g., SUSY) models give predictions for a p (and a n ) which can be next directly 
compared with experimental results. 

In the SI case the conversion is straightforward. This is because now the enhancement 
factor is proportional to the square of the atomic number, Ca = —ttt [Z f p + (A — Z)f n ] 2 

where f p and f n are the effective WIMP couplings to protons and neutrons, respectively. For 
Majorana WIMPs f p ~ f n and so one typically has Ca/C p ~ Ca/Cu ~ A 2 and the conversion 
does not depend on the specific WIMP type. For massive Dirac neutrino-like WIMPs f p ~ 
and Ca/C p ~ Ca/C u ~ (A — Z) 2 and again the conversion does not depend on the WIMP 
model or its parameters. 

In the SD case however the situation is more complex. Here the enhancement factoiQ is 
given by 

8 J + 1 

C A = - (a p (S p ) + a n (S n )) 2 — — , (3) 

7T J 

where a p and a n are (WIMP-type dependent) effective WIMP-proton and WIMP-neutron 
couplings and (S Ptn ) = (N\S Ptn \N) are the expectation values of the proton and neutron spins 
within the nucleus and J is the total nuclear spin. In the particular case of free nucleons one 
finds C Pin = ^ a 2 n and a Pi7l = ^■G 2 F ^L 2 n a^ ) n where (S p ) = 0.5 = (S n ) has been used. 

1 Note that the enhancement factor Ca as defined here (and in Ref. ||) differs slightly from a similar 
quantity I a used in Refs. 14 0|: Ca = -Ia- 
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This definition of Ca normalises the spin-dependent nuclear form- factor F 2 (q) used in 
calculating nuclear recoil energy spectra to unity at q = 0: 



where 



5(g) = c^5 o(g) + a\Sn{q) + a aiS i(q) (5) 
= (a P + a n ) 2 S 00 (q) + {a p - a n ) 2 S n (q) + (of,- a^j S i(q). (6) 

Here the 5jj are the respective isoscalar, isovector and interference term form-factors for 
nucleus N (assumed to be known from nuclear calculations) and oo = a p + a n and a\ = 
dp (in are isoscalar and isovector coefficients. Using this definition most of the WIMP 
model-dependencies in S(q) are absorbed into Ca, leaving F 2 (q) and hence the shape of recoil 
energy spectra relatively model- independent. Recent calculations || have shown that F 2 {q) 
still contains some residual model-dependencies owing to differences in the g-dependence of 
the Sij form-factors. These differences are however small for most nuclei of interest, including 
Na, I and F ||, and will henceforth be neglected. 

It is clear from Eqn. (^) that, in the SD case, converting WIMP-target cross section limits 
a\ m to the WIMP-proton cross section Op becomes problematic. This is because the 
enhancement factor Ca now receives contributions from both proton and neutron terms. As 
we will argue below, due to the presence of WIMP-dependent coefficients a Pj „, both of these 
contributions can be of comparable order and even of opposite sign. As a result, the ratio 
Ca/C p in Eqn. @, and hence the derived limit <7p im ^ will depend on the assumed type of 
WIMP. 

In the early (single-particle or odd-group model) calculations the nuclear spin was always 
assumed to be dominated either by the proton or by the neutron term in Eqn. (^). Thus 
when dealing with odd-proton targets, such as Na or I, the a 2 in the expression for Ca was 
conveniently cancelled by the a 2 in the analogous expression for the WIMP-proton cross section 
enhancement factor C p . Model dependencies were thus eliminated. 

The conversion is however complicated when dealing with odd-neutron targets because 
of the factor (a p /a n ) 2 remaining in the expression for ap 1 ™^ in Eqn. @. This ratio is not 
guaranteed to be constant and in general is WIMP model-dependent. In the important case of 
SUSY neutralino WIMPs early estimates of the Aq values used in calculating a p and a n were 
nevertheless such that the ratio <J p /a n = (a p /a n ) 2 was always of order unity, independent of 
the neutralino composition. Later estimates of Aq [Q] have however shown that although this 
is still true for predominantly higgsino neutralinos, the ratio for gaugino neutralinos is highly 
SUSY model-dependent and can vary by several orders of magnitude (as demonstrated by 
Fig. H for models from the database described in Refs. 0, ^)- Gaugino- like neutralinos as DM 
WIMPs are strongly favoured by a combination of naturalness and cosmological arguments || 
in which case this problem becomes particularly acute. 

A still further problem arises when using more recent shell-model calculations for (S p ) 
and (5 n ) (H, HJ. These indicate non-zero contributions to the nuclear spin from both protons 
and neutrons, and in this case the a 2 factor cannot be cancelled from even odd-proton nuclei. 
Hence even if one of these contributions is larger than the other, as is often the case, then the 
WIMP-dependent ratio of a p and a n can be such that both contributions to the cross section 
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are substantial. Furthermore, a p (S p ) and a n (S n ) can in general be of opposite sign and similar 
magnitude. Hence, a more proper way of writing the enhancement factor would be 

c A = - (K(S p )\ ± |a n (S n )|) 2 ^±1. (7) 

TT J 

In general therefore, depending on the particular WIMP type, the WIMP-target cross section 
can be considerably reduced relative to that of its constituent nucleons. In these circumstances 
limits set by assuming the simple case of higgsino neutralino or heavy neutrino WIMPs (con- 
stant (a p /a n )) would prove to be unduly optimistic. 



3 An Alternative Procedure 

We will now present a new method for deriving limits on spin-dependent WIMP-nucleon 
interactions. The method will be free from the problems described above and will allow 
one to derive experimental limits in a way which is independent of the type of the assumed 
(Major ana) WIMP. 

We will start by identifying the separate proton and neutron contributions to the total 
enhancement factor Ca 

C\=*{a p (S p )f J -±±, Cl=*{a n {S n )fl±± (8) 

TT J TT J 

In light of Eqn. @, this gives C A = (yCf ± V^a) • Following Eqns. ® and (|) we thus 
define the proton and neutron contributions a* and a\ to the total cross section a a as: 

a\ = AGWaCa < = A&fAPI- (9) 
Using Eqns. (|l]), M) and (g), one can express a a as 

a A =(M±\^X) ■ (10) 



We note that o~ A and a\ are not measured cross sections. They are nevertheless convenient 
auxiliary quantities which identify separate proton and neutron contributions to the total cross 
section oa- 

We will now proceed as follows. (See the Appendix for a more rigorous treatment.) We will 
first make an auxiliary assumption that a a — o~ A . In other words, we assume that the total 
WIMP-nucleus cross section is dominated by the proton contribution only. We then define the 
WIMP-proton cross section limit a p m ^ corresponding to the WIMP-target A cross section 
limit a\ m as 

u 2 1 

lim(A) lim x 

Analogously, the WIMP-neutron cross section limit cr^ 111 ^^ is defined by assuming that aA — 
a\ and writing 



^^/piF' (12) 

rA ^Al^n 
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It is clear that the use of the ratios C P A /C p = 4/3(S p ) 2 (J+l)/J and C n A jC n = 4/3(S n ) 2 (J+ 
1)/J ensures the cancellation of the WIMP-dependent a 2 and a 2 terms contained within the 
WIMP-target cross section a a and hence ensures WIMP model-independence. Values of 
C P A /C P and C A /C n for typical nuclei of interest obtained using data from Refs. g |. || are 
listed in Table g. 

These properties can now be used in expressing WIMP-independent experimental limits 
on the WIMP-nucleus SD interaction cross section a\ m in terms of a p and a n which are the 
quantities whose values are predicted by specific theoretical models. If an experiment publishes 
the limits a p m ^ and then Eqn. ([n]) can be used to define a WIMP-independent 

excluded region in the o p - a n plane as 




> 1. (13) 

Because of relative sign ambiguity, the condition (^) implies two bounds corresponding 
to constructive and destructive interference. A conservative limit corresponds to the relative 
minus sign which reduces the overall WIMP-target cross section. We note, however, that in 
comparing with a specific theoretical (e.g. SUSY) model there will be no sign ambiguity: the 
theoretical model predicts not only a p and a n but also the signs of a p and a n . In this case the 
sign in Eqn. ( |T3| ) is known and is given by the sign of (a p (S p ))/(a n {S n )). 

An alternative way of expressing the limits in this procedure is to consider exclusion regions 
directly in terms of the fundamental WIMP-nucleon coupling coefficients a p and a n . In this 
case Eqn. (O) is replaced by 




24G> 2 ' (14) 

where we applied Eqn. (Q) to the case of the nucleons to obtain a p = 24Gpfi p a p /TT and 
a n = 24G 2 7 /i 2 a 2 /7r. The small proton-neutron mass difference has been ignored. 

The variables a p and a n can have either sign, and the relative sign inside the square is 



now determined by the sign of (S p )/(S n ) only. Eqn. (14) corresponds geometrically in the a p 
- a n plane to excluding a region exterior to two parallel lines whose slope has opposite sign to 
{S p )/(S n ). There is no limit on a p or a n between these lines. This region extends to infinity 
in both directions. 

So far, we have presented the method for one species of target nucleus only. A gener- 
alisation to two or more nuclei in the same target is straightforward. Analogously to the 
WIMP-proton cross section limit in the current method 0, the limits a p m ^ A ^ and o^^ 1 ' 
from different nuclei Ai in the target A can be combined by using 

lim(A) = X! lim(Ai) ' lim(A) = X! WA~) ' 

O p Ai Cp Cn Ai &n 

It should be noted that when calculating the combined excluded region in a p - <r n parameter 
space for two or more nuclei in the same target material it would be incorrect to use Eqn. (|i"3 ) 
with the combined limits <rp im ^ and <r^ m ^ calculated using Eqn. ([11]). The correct approach 
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is to use instead the generalisation of Eqn. ([b^) given by 




0V, 



Um(Ai) \/ lim(Ai 
Up | On 



> 1. (16) 



A similar generalisation can also be applied to Eqn. (|l~4|). In this case for given cr ; 



lim(A) 



V 



and an optimum limits on the coupling coefficients a p and a n could be obtained by using 
two different target nuclei with (S p )/(S n ) of opposite sign, such as are found in NaCl or NaF. 
The allowed region would then lie inside the intersection of the two bands of opposite slope. 

As a practical illustration, we consider limits from a synthetic data set (assumed to be 
from a Nal detector) using both the current and the proposed techniques. The data set 
consists of recoil energy dependent event rate limits which have been converted to WIMP 
mass dependent target nucleus cross section limits using nuclear kinematics and Eqn. (Q) as 
described in Ref. ||. 

In Fig. |2](a) we plot the limits calculated using the current technique (i.e. using the ex- 
pression for <7p im ^ given in Eqn. (||) with Ca from Eqn. (||)). We do this for three different 
neutralino WIMP compositions (i.e. different ratios a p /a n from Fig. ||). This results in over- 
all (combined Na and I) limits <7p im ^ NaI ^ which, for a given WIMP mass but different type, 
can be different by almost two orders of magnitude : a thick solid curve corresponds to the 
limit for the higgsino WIMP (a p /a n ~ 1.5), which is currently commonly assumed, while the 
two thick dash-dotted lines correspond to the two gaugino WIMP cases giving the extremal 
limits for WIMPs of mass 100 GeV/c 2 . These two curves illustrate the considerable effects 
of destructive (upper curve) and constructive (lower curve) interference between proton and 
neutron contributions to the nuclear spin. 

Results obtained using the proposed technique are presented in Figs. |2|(b) and (c). We 
plot the limits <7p im ^ and a 1 ™^ calculated using Eqns. (11) - (|l~2|) and hence avoiding the 



assumption of a specific WIMP composition. We can clearly see that the limits in windows 
(a), (b) and (c) differ considerably. In all the cases at small WIMP masses the best limits are 
provided by Na due to its slowly- varying form- factor F 2 (q). However, at larger masses iodine 
provides a better limit due to its larger value for C A / C p (Table [j]) . The spin of both nuclei 

is carried predominantly by protons and so for this target the limits <7p im ^ are in general 

lim(A) 
superior to a n 

The individual limits o~ p im ^ and On presented in Figs. |2| (b) and (c) contain all the 
information which allows one to draw, for a given WIMP mass, exclusion regions in the a p 
- a n plane by using Eqn. (|i"6|). This is shown in Fig. |J| We use the same data set as before 
and set a WIMP mass of 100 GeV/c 2 . Figs. |3|(a) and (b) correspond to minus and plus signs 
respectively in Eqn. fll~6|), although in the absence of knowledge of this sign the former plot 
gives the more conservative limits. 

Finally we note that, in the presented technique, one can incorporate limits a^f^' set 
on the spin-independent WIMP-nucleon cross section o~si by constituent nuclei using the 
technique of Sec. 1: 



E 




> 1- (17) 
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In summary, in contrast to the currently used practice of using only the WIMP-proton 
cross section a p for presenting the SD WIMP-nucleus cross section limits, this alternative 
method effectively uses both the WIMP-proton and WIMP-neutron cross sections a p and a n . 
This ensures that the WIMP-model dependence of the experimental SD cross section limits is 
removed. The method applies to targets with one or more species of nuclei in the target and 
allows unambiguous comparisons with theoretical predictions. 

4 Discussion and Conclusions 

In contrast to the currently used procedure, the new technique makes it possible for a given 
experiment to set limits on the spin-dependent WIMP-nucleon cross section in a WIMP- 
independent way. Another significant advantage of the new technique is that it makes the 
direct comparison of experimental results with theoretical predictions possible. Given the 
WIMP-proton and WIMP-neutron exclusion plots from a particular experiment, one can de- 
termine whether a specific choice of parameters is allowed or excluded for a given WIMP mass 
simply by using Eqn. (|l~7|). In fact, all the information needed to allow or exclude any WIMP 
candidate (like SUSY neutralino WIMPs with arbitrary composition, non-SUSY heavy neu- 
trino, etc.) can be presented in only three figures (e.g. Fig. §(b) and Fig. |2|(c) plus one SI 
cross section limit plot). No additional plots similar to Fig. || need be drawn to determine the 
exclusion region for a particular WIMP mass. 

It is not surprising that dark matter experiments using odd-proton targets will generally 
give very good WIMP-proton cross section limits and rather weak WIMP-neutron cross section 
limits for the same WIMP-target cross section limits. Just the opposite will be true with 
experiments using odd-neutron targets. This is what one should in fact expect: odd-proton 
and odd-neutron targets are effectively measuring two quite different and unrelated quantities 
(a p and a n ). It is only in the particular case of higgsino neutralino and heavy neutrino 
WIMPs H| that the limits from such targets calculated using the new technique are still 
comparable, by using a constant value for a p /a n ~ 1.5 to combine Op m ^ and On into a 
single limit on the WIMP-proton cross section. These limits are equivalent to those obtained 
using the current technique (i.e. Fig. ||(a)). In any case, experiments may find it useful to 
publish such limits in addition to the a p and a n limit curves calculated using the new technique 
(i.e. Fig. |(b) and Fig. |(c)). 

In the event of a discovery of spin-dependent WIMP-nucleon interactions it is likely that 
a different procedure from that described above would be required to analyse the data. In 
this case a fit would likely be performed to the observed nuclear recoil energy spectrum, with 
parameters such as the WIMP mass and cross sections on each target nucleus being considered 
to be free. This would in turn define an allowed region in & p - & n parameter space for each 
allowed WIMP mass and it would be this region which could be compared with theoretical 
predictions. 

In conclusion, a new technique has been presented for deriving spin-dependent WIMP- 
nucleon cross section limits from direct detection dark matter experiments. This technique 
retains the attractive features of the current procedure, including the ability to combine limits 
from individual target nuclei. The new technique, however, has the additional advantage 
of allowing the calculation of WIMP-nucleon spin dependent cross section limits in a WIMP- 
independent manner and of making it possible to compare experimental results with theoretical 
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predictions. 
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Appendix 



Here we present the basic formalism behind the alternative method presented in Sec. 3. For 
the sake of completeness, we will include here some expressions given already there. 

The separate proton and neutron contributions to the total enhancement factor Ca given 
in Eqn. (||) are given in Eqns. (||) 



C p A = -(a P {S p )f±±±. 

7T J 



C: 



- (a n {s n )y 



7T 



J + l 
J 



. Following Eqns. (|j) and (||) the proton and 



so that, from Eqn. ([?]), Ca 



neutron contributions a v A and a A to the total cross section a a are defined as (Eqns. (||)): 



Using Eqns. ([[]), (0) and @, one can then write a a 



(Eqn. ©). 



Applying Eqn. ([l]) to protons and neutrons respectively and comparing with Eqns. (||) we 
find the following relations 



o,, 



Va <-m 



a 2 r 

,,2 rm ' 



(18) 



It is clear that the use of the ratios C P /C A and C n /C A ensure the cancellation of the ai and 
terms contained within WIMP-target cross section oa and hence ensure WIMP model- 
independence. Unfortunately, the quantities a A and a A cannot be measured directly. We will 
therefore make independently the assumptions that 



OA 



cr 



A' 



OA 



in which case Eqns. become 



a 2 c p 
Va °a 



OA ',2 nP ' 



o„ 



OA 2 



2 nn ■ 
I 1 A A 



(19) 



(20) 



In other words and are the WIMP-proton and WIMP-neutron cross sections derived 
from the WIMP-target cross section cja by assuming that it is dominated by its proton and 
neutron contributions respectively. These new quantities are related to a p (a n ) and a A (o A ) 
by 

(21) 



Op 

p 



cr 



A 
OA 



On 
A 



O, 



OA 



Hence from Eqn. (|To|), 




^ =1- 



(22) 



Experiments set limits a A m on a a, while a p and a n are quantities whose values are predicted 
by theoretical models. It is our goal to set limits on a p and a n by using the measured values 
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of c^ m . In direct analogy with Eqn. @ we then define the WIMP-proton cross section limit 
^hm(A) correS p 0nc jj n g to the WIMP-target A cross section limit cr^ m as 

lim(A) _ lim Mp &P_ rr lim(A) _ „lim Mn fit 

A 4 a °a Ma °a 



which are Eqns. (O) - (12) in Sec. 3. If an experiment publishes the total SD WIMP-target 



cross section limit a^ 11 in terms of the corresponding limits <7p im ^ and Cn™^ then a WIMP- 



independent excluded region in the a p - a n plane will be given by the condition ( |L3] 
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Tables 



Table 1: Values of (S p ), (S n ), C P A /C P and C\jC n for various nuclei. Values of (S p ) and (S n ) 
for Na, Te, I and Xe are taken from Ref. H|. Values for F are taken from Ref. ||. All others 
are from the review of Ref. § and the references contained therein. 



Figures 

Figure 1: The ratio c p /a n plotted against neutralino composition Z g /(1 — Z g ) for models 
from the database built up in Refs. |, g. Here Z g is the gaugino fraction and (1 — Z g ) 
is the higgsino fraction. Small values of Z g /(1 — Z g ) correspond to predominantly higgsino 
neutralinos while large values correspond to predominantly gaugino neutralino. In all models 
plotted, the neutralino is a good dark matter candidate (i.e. its relic density is in the range 
0.025 < Q x h 2 < 1). The neutralino-proton and neutralino-neutron cross sections o p and a n 
are calculated as in Ref. j7|]. 



Figure 2: Exclusion regions for simulated data from a Nal detector. Figure (a) shows in 
the current technique the limits <7p im ^ calculated using Eqn. ( |l5|) for the case of three differ- 
ent neutralino WIMPs. The thick solid curve is the combined limit for a higgsino-like WIMP 
(corresponding to a p /a n ~ 1.5), which is currently commonly assumed. The two thick dash- 
dotted curves are the combined limits calculated for two different gaugino neutralino cases 
assuming either destructive (upper curve) or constructive (lower curve) interference respec- 
tively between proton and neutron contributions to the nuclear spin. Also shown are the 
individual limits from Na (dashed) and I (dotted) nuclei contributing to the combined limit 
on the higgsino-like WIMP. 

Figures (b) and (c) show the limits <7p im ^ and respectively calculated from the same 

data using Eqns. (JTl|) - ( Jl2| ) in the framework of the new technique. Dashed (dotted) curves 
again correspond to individual limits from Na (I) nuclei and thick solid curves to combined 



limits obtained using Eqn. 15. The limits are neutralino WIMP type independent. 



Figure 3: Exclusion regions in a p - a n plane are plotted in the case of 100 GeV/c 2 WIMPs 
calculated from the data of Fig. |2|(b) and (c) using Eqn. [H]. Figure (a) is for the conservative 
case of destructive interference ((a p (S p )) / (a n {S n )) < 0) while Figure (b) is for the case of 
constructive interference ((a p (S p ))/(a n (S n )) > 0). Dashed (dotted) curves correspond to 
limits from Na (I) alone while the full thick curves show the combined limits. Note that 
in Figure (a) no limit can be set by any one nucleus when cr p / = o" n /<7n m ^ due to 
destructive interference between the proton and neutron contributions. Combination of limits 
from two nuclei with different <7p im ^ /a™^ allows such a limit to be set however. 
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